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Abstract 

We consider a network model, embedded on the Manhattan lattice, of a quantum localisation 
problem belonging to symmetry class C. This arises in the context of quasiparticle dynamics 
in disordered spin-singlet superconductors which are invariant under spin rotations but not 
under time reversal. A mapping exists between problems belonging to this symmetry class 
and certain classical random walks which are self-avoiding and have attractive interactions; we 
exploit this equivalence, using a study of the classical random walks to gain information about 
the corresponding quantum problem. In a field-theoretic approach, we show that the interactions 
may flow to one of two possible strong coupling regimes separated by a transition: however, using 
Monte Carlo simulations we show that the walks are in fact always compact two-dimensional 
objects with a well-defined one-dimensional surface, indicating that the corresponding quantum 
system is localised. 

Short Title: Localisation and the Manhattan lattice 

PACS numbers: 72.15.Rn, 05.40.Fb, 64.60.Ak, 05.50.-hq 
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1 Introduction 



Disordered electronic systems, and the associated metal-insulator transitions, have generated much 
interest in recent years. In this paper, we consider a two-dimensional system of non-interacting 
particles which is invariant under spin rotations but not under time reversal. According to the 
classification set out by Altland and Zirnbauer such a system belongs to symmetry class C. This 
class is one of several which are distinct from the three Wigner-Dyson classes and have realisations as 
Bogoliubov-de Gennes Hamiltonians for quasiparticles in disordered superconductors. In particular, 
class C arises in spin-singlet superconductors with broken time-reversal symmetry for orbital motion 
but negligible Zeeman splitting 

For two-dimensional systems with broken time reversal symmetry, whether from the Wigner- 
Dyson unitary class or from class C, there exists the possibility of a delocalisation transition of 
the quantum Hall type. Recent studies of network models for systems belonging to class C 01 
1^ m have focussed on this transition, yielding several interesting results. Notably, one such 
network model describes the spin quantum Hall effect, with two insulating phases separated by a 
delocalisation transition analogous to the quantum Hall plateau transition. The two phases exhibit 
quantised Hall conductance differing by an integer. It has been shown [7j that this quantum 
problem may be mapped to that of the hulls of classical percolation clusters in two dimensions, a 
well-understood problem. The particular network model used in these studies was based on the 
L-lattice (see Fig. ^), and is a natural generalisation of the Chalker-Coddington modeljH] which 
was designed to study the ordinary integer quantum Hall transition. 

The L-lattice is a square lattice consisting of directed edges and nodes, such that every path 
following the directed edges must turn left or right at each node. It possesses a particular symmetry 
under inversion about a diagonal axis which reflects the symmetry of the corresponding quantum 
hamiltonian about a particular energy. 

However, in [Jj it was shown that the mapping of network models in class C to classical random 
walk models is more general, and holds on any directed lattice (at least if each node is of degree 
4.) These walks are in general trails, in that they may not pass over a given link more than once, 
but they also have attractive interactions at each node. In the mapping, localised and extended 
phases of the quantum system map to regions of parameter space where the corresponding classical 
walks either close after a finite number of steps (almost surely) or they escape to infinity. The 
localisation behaviour exhibited by a given network model, and the nature of the corresponding 
classical problem is expected to be strongly dependent on the structure of the underlying lattice. 
Therefore, by studying classical problems formulated on different lattices, we might uncover dif- 
ferent types of localisation behaviour, resulting in the competition between the edge exclusion and 
nodal attraction. 

The classical random walks in which we are interested, together with other kinetic self- avoiding 
trails, correspond to the classical scattering of light by random arrays of mirrors laid out on a square 
lattice. Problems of this type have been studied extensively [SlEl^nEl. These processes may 
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(a) the L-lattice 



(b) the Manhattan lattice 



Figure 1: The L and Manhattan lattices 



also be viewed as history-dependent kinetic random walks. For example, in terms of the mirror 
model, the walker lays down a mirror at random on visiting a node; on return to that node, it 
finds the mirror already fixed there. Therefore, models of this type may be classified as "true" 
self-avoiding walks. This class of walks has been studied in the limit of weak scattering using 
field-theoretic renormalisation group methods |141 1151 [T^ ; we shall make use of such techniques in 
section 121 of this paper to study such walks. We note that, as for Anderson localisation, the critical 
dimension for these processes is two. 

As mentioned above, the network model describing the spin quantum Hall effect is formulated 
on the L lattice (Fig. ^). In the heuristic justification for studying this lattice in this context (and 
for the original Chalker-Coddington model) it is seen as a distorted version of a random lattice, 
whose plaquettes correspond to the regions where locally either — y is < or > 0, where E is the 
Fermi energy and V is some slowly varying random potential. In the semi-classical approximation, 
there are edge states propagating along the curves which divide these regions. Along these, the 
quasi-particle wave function suffers a constant spin (or phase) rotation. On distorting the random 
geometry to the L-lattice, these become quenched random link variables of the model. 

In this paper, we consider a different lattice, the Manhattan lattice. It shares with the L-lattice 
the property that at each node there are 2 incoming and 2 outgoing edges or links, so that the 
theorems of apply. On this lattice, a walker may either continue straight on at a given node with 
a given probability or make a turn. The allowed directions for turning at each node alternate right 
and left, as shown in Fig. The overall lattice is therefore regular, with a similar size unit cell to 
the L-lattice. In the network model, at each node there is a real unitary S'-matrix which describes 
scattering from the incoming edges to the outgoing ones. Just as for the models on the L-lattice, 
this is in principle also a quenched random matrix, but for simplicity we shall take it to be the 
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same at all nodes. Thus the 5-matrix elements for continuing straight on at each node are taken 
to be cos 9, and those for turning are it sin 0. When cos 6 = 1, the quasiparticles move ballistically 
along the lattice axes, while for sin 9 = ±1 they move around the perimeter of single plaquettes. 
The existence of a ballistic limit distinguishes this lattice from the L-lattice. On the other hand, 
the Manhattan lattice appears to enjoy no special symmetry, as does the L-lattice about = 7r/4. 

The general results of show that the single-particle Green's function and the two-point 
conductance for a class C model on this lattice are related to the properties of classical trails, or 
random walks which visit each edge no more than once, and, on first visiting a given node, turn 
L or R with probability p = sin^ 9, or continue straight with probability 1 — p = cos^ 9. For the 
L-lattice, the limits p ^ and p — > 1 show clearly two distinct behavioural phases, characterised 
by different edge states; a delocalisation transition occurs at p = 1/2. For the model formulated on 
the Manhattan lattice, however, we might expect to find different behaviour, since it is clear that 
there are no corresponding edge states as p — > 0, 1. We note also that a kinetic self-avoiding trail 
on the Manhattan lattice may cross itself with probability p, but this is never permitted on the L 
lattice. 

At p = 1, the walks are simply loops defining the elementary plaquettes of the lattice. In fact, 
it may be argued jlfij that the trajectories are loops of finite extent at least for p > ^: if we place 
a mirror at 45° at each of the vertices where the path turns through ±90°, they may be thought 
of as forming clusters of bonds on a square lattice at 45° to the original. For p > ^ the mirrors 
percolate: that is, they form an infinite cluster which crosses the entire lattice, with probability 
one. The trajectories are confined to the voids in this cluster, which are finite. For < p < ^, 
this argument is inconclusive: the trajectories will scatter from the finite clusters of mirrors, but 
this may or may not lead to classical localisation. At p = all trajectories are straight lines; the 
states to which they correspond are spatially extended. Therefore either almost all trajectories are 
localised for all p > 0, or there is a transition at some finite value < Pc ^ from extended to 
localised behaviour. 

In this paper, we use two approaches to study the kinetic processes described above. In section 
m we develop a field-theoretic description of the weak-coupling regime of this particular type of 
"true" self-avoiding walk. We begin by estimating the diffusion constant for the Manhattan lattice, 
and show that the problem is in this regime as p ^ 0. We then use field-theoretic renormalisation 
group methods to analyse the general class of such problems: the RG trajectories always run away 
from the weak-coupling fixed point, with possibility of two distinct strong-coupling behaviours 
separated by a critical trajectory. However, in section |31 we use Monte-Carlo methods directly to 
simulate kinetic self-avoiding trails on the Manhattan lattice. These simulations are only practical 
for values of p larger than about 0.2. However, they all indicate that the walks are localised, 
with a mean length and spatial extent which decrease with increasing p. These numerical results 
thus complement those of section |2j In addition to their relevance to the corresponding quantum 
localisation problem, these classical processes are of interest in their own right. In section 01 we 
discuss possibilities for related future work. 
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2 Analytic results 



2.1 Estimation of the diffusion constant 

As a first step in analysing the weak-interaction limit of this model, we estimate the diffusion 
constant of a random walker on the Manhattan lattice, subject to no restrictions on return to the 
same node or link. A diffusive path r(t) which takes a step of finite length at each unit time exhibits 
the behaviour 

{[rit) - rm') = 4Dt , (2.1) 

where the lattice-dependent quantity D is the diffusion constant. For convenience, we shall take 
r(0) = in what follows. The diffusion constant for a given random walk may be calculated by 
recasting the position coordinates of points visited on the path as a sum rt of complex numbers: 

rt = zi + Z2 + ... + zt, (2.2) 

where Zj^i = Zj e*^^ . For the random walks in which we are interested, we require additionally 
\zj\ = a for all j, where a is the lattice bond length. The distribution of 9j depends on the details 
of the lattice; the 9j are assumed to be uncorrelated, so that each step of the random process 
is independent of all previous history. For the Manhattan lattice, the distribution for the 9j is 
somewhat complicated because the allowed turning directions at a given node are dependent upon 
the position of the node on the lattice. For turning probability p, the mean free path of a particle 
moving randomly on the Manhattan lattice Tj-Qfp ^ 1/p; hence for small p, r^fp » 1. We therefore 
assume that the particle is unable, on the length scale of the mean free path, to tell whether it will 
next change direction at a left-turn or right-turn node. Thus the options of turning left or right 
are assigned equal probability in this approximation. The corresponding distribution for 9j is 

with probability 1 — p 
-7r/2 with probability p/2 (2-3) 
+7r/2 with probability p/2 , 

with angles measured in an anti-clockwise direction. Thus (e*^^ ) = 1 — p + e*'^/^ p/2 + e^*'^/^ p/2 = 
1 - p for all j, and (e'(^i+-+^'-)) = (e'^i)'' = {\ - vY . Now from ^J^, we have 

n = (1 + e^^i + . . . + e*(ei+-+e*-i)) . (2.4) 

Hence, since (|r(t)p) = (nr*), 

(|r(i)|2) = {i2-p)p-h - 2(1 -p)p-2 (1 - (I-pY)) . (2.5) 
At large t, {I — pY —>■ 0, yielding 

(|r(t)|2)p.|zi|2((2-p)p-it-2(l-p)p-2) . (2.6) 
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Thus for the Manhattan lattice, 

D^{2-p)p-K (2.7) 

We note that at p = 0, the diffusion constant becomes infinite; this characterises superdiffusive 
behaviour. In this particular example of superdiffusion, we also find that the mean square dis- 
placement {[r{t) — r(0)]^) diverges faster than linearly with time. In fact, it is easy to see that 
{[r{t) — r(0)]^) = since for p = the random walk proceeds in a straight line. We therefore 
see a transition to a different kind of behaviour exactly at p = 0. 



2.2 A field theoretic approach 

2.2.1 Peliti's field theory for true self-avoiding walks 

The starting point for our renormalisation group approach is the field theoretic description proposed 
by Peliti |14| II 5| I13j for the class of true self-avoiding walks to which our bond-avoiding paths on 
the Manhattan lattice belong. The theory is composed of two scalar fields ip{r, t) and ip{r, t), which 
respectively create and destroy walks. The Hamiltonian Ti is given hy 7i = TCq + 7i[, where 7io is 
the 'free' Hamiltonian 

d 



fOO f 

no = J dt J d'^ri!{r,t) 



,p(r,f) (2,J 



and Til is an interaction Hamiltonian whose precise form depends on the nature of the process 
under consideration. D represents the bare diffusion constant, as computed above. From this 
class of Hamiltonians arise well-defined Feynman rules for the theory. The correlation function 
G{r,t) = (■0(0, 0) V'Cr, i)) yields the probability density for a walk beginning at the origin to be 
found at position r after time t. After applying a Fourier transform to the spatial components of 
G{r, t) and a Laplace transform to the temporal component, we obtain the bare propagator 

G'(p,/x) = Go(p,/x)^(ai + I?/)-^ (2.9) 
The generalised interaction contribution is given by 

7(p, q, Pi) = (pi-q) 91 + (pi-(pi + q)) 92 + (pi-p) 93 , (2.10) 

where p is the earliest incoming wavenumber and pi is the latest outgoing wavenumber. The 
values of the couplings gi depend on the nature of the random walk. For example, for Peliti's true 
self-avoiding walk, one should take 91 = g and 92 = 93 = 0; this corresponds to the interaction 
Hamiltonian 



ni = -g dtdt' J d'^r (^^/;(r, t')Vi^{r, t') 



, (2.11) 

Since the coupling constants are dimensionless quantities, the upper critical dimension of this class 
of theories is dc = 2. The general theory may be renormalised by dimensional regularisation and 
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minimal subtraction in two dimensions. Mass renormalisation is trivial. However, the diffusion 
constant is renormalised by a factor Z so that D = Z Dr, and renormalised coupling constants are 
given by Ui k^, i = 1, 2, 3, where k is the renormalisation wavenumber. A calculation to order one 
loop gives diffusion constant renormalisation factor 

Z = 1 + e-\u3 - ui) (2.12) 



-j^ = —eui + |uf + — 2uiU2 — SuiU3 (2.13) 



and renormalisation group equations 

du 
"d^ 

— = —€U2 — uf — U2 + 3uiU2 + U1U3 — 3n2n3 (2-14) 
dr 

^7 = —eu3 — + 2uiU3 — U2U3 (2.15) 
where r = log(«;'/«;) is the logarithm of the scale parameter. 

2.2.2 Constructing an interaction Hamiltonian 

In order to apply these general results to the random walks in which we are interested, we now 
examine in detail the types of interaction which might arise for kinetic self-avoiding processes on 
the Manhattan lattice. There are two competing interactions involved: the node interactions are 
attractive, and clearly the bond interactions in this type of self-avoiding process must be repulsive. 
Thus one way to look at the effect of these interactions is to begin with a free kinetic walk on the 
Manhattan lattice and then introduce appropriate weights to model self-avoiding characteristics. 
We impose a repulsive weight exp(— /?FM(M— l)/2) (with F > 0) for every bond visited M times 
and an attractive weight e-xp{+pV'N{N — l)/2) (with V' > 0) for every node visited N times, and 
examine the resultant modification to the free kinetic walk. To work to first order in V and V, it 
is sufficient to consider only cases where a path has visited each node and bond on the lattice no 
more than twice. 

On average, there is no directional bias to these processes on the Manhattan lattice; therefore, 
no drift terms, or equivalently no first order gradient terms, should appear in the interaction 
Hamiltonian. However, higher order gradients may be obtained by examining the hopping terms 
generated as the particle moves on the lattice. These should respect the global rotational invariance 
of the kinetic walks. Consider, for example, a move from ri to r2 at time t2- This is effected by the 
probability-conserving term 'tp{r2) ip{ri) — tp^ri )'0(ri). Note that the particle density at position 
r and time t is given by 'tp{r, t) 'tp{r, t), and takes values or 1. Now suppose that the midpoints of 
the exit bonds at a particular node have position vectors ri, with i = 1, . . . ,4. Thus to describe a 
first visit to the node at time ti and a subsequent visit at time t2, approaching from ri and leaving 
via r2, we postulate an interaction term of the form 

[V'(r2, t2) '0(ri, t2 )-'!/' (ri, t2) V'(ri, i2)] 

4 (2.16) 



7 



1 




1 


< 

1 


< 

3 4 

(x,v) 




— 

/ 


5» 

\ ^ 





Figure 2: A unit cell on the Manhattan lattice with centre at the point {x,y); the possible node 
types are labelled 1 to 4. Each lattice bond has length 2a. 



with the ai appropriately chosen to give the correct Boltzmann weights. Suppose that the possible 
routes through a given node are ri ^ r2, ri ^ rs, r4 ^ r2 and r4 ^ rs. Then applying weights 
V and V as defined above, if the path passes from ri to r2 on both visits to the node, we obtain 

ai + a2 = -2V + V' . (2.17) 

Considering each of the other possible routes through the node, we find also 

ai + a3 = -F + F' (2.18) 
a2 + a4 = -V + V' (2.19) 
03 + a4 = +V' (2.20) 

We therefore have a set of simultaneous equations for the Note that these are invariant under 
the transformation 01^4 — > 01^4 + A, 02,3 — > 02,3 — ^ for any A, so that there is always some 
arbitrariness in the solutions. However, if V' = 0, i.e. there is no node interaction, then ai = 
a2 = —V, a3 = 04 = is a solution. Similarly, ii V = 0, i.e. there is no bond interaction, then 
«! = 02 = 03 = 04 = V' /2 is a solution. We thus obtain bond interactions of the form 



-V (^(rj, i2) V'(ri, t2) - ^(ri, t2) V'(ri, ^2)) 

X (Vi(rj,ti)^(ri,fi) + 7/;(rj,ti)V'(rj,ti)) 



(2.21) 



and node interactions of the form 

1 

2 



-V (^(rj, f2) V(ri, t2)-^(ri, t2) V'(ri, t2)) 



^(r',ti)V(r',ti), 

r'Snode 



(2.22) 



where, on the second visit, the path approaches the node from rj and leaves via rj. 

A unit cell on the Manhattan lattice contains four different node configurations, as illustrated in 
Fig. [21 In order to construct the interaction terms in the Hamiltonian then, we must compute the 
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contribution for both node and bond interactions from each of the nodes in the unit ceh. Consider, 
for example, node 1 in Fig. |21 positioned at (x — a,y — a). A possible route through this node 
might arrive at the point = (x — 2a, y — a) at time t2 and leave via rj = (x,y — a). The second 
order gradient expansion for this hopping process is given by 

-0(rj, t2)'0(ri, ^2) - i'iri, t2)V'(ri, ^2) 
= i){x, y- a, t2)'>p{x - 2a, y, 12) 

-Tpix -2a,y,t2)i'{x -2a,y,t2) 

(2.23) 

= 2a ^ V2 - 2a2-^V2 - 2a^-^ij2 
ox ax^ oxoy 

dx dx dx dy 

In the last line above, the subscript 2 indicates time t2 and the position argument (x, y) of the ip, ip 
has been omitted for clarity. Combining the contributions for all the possible routes through each 
of the four nodes belonging to the unit cell, we obtain the node interaction term 

-16^a2(3(VV;(r,t2).V^(r,t2))(V5(r,ti)V'(r,ti)) ^ 

2 (2.24) 

+ 4V.(a,(VV;(r, t2))V'(r, t2)^(r, ti)V(r, ti))) 

and bond interaction term 

8Fa2(3(V^(r,t2).VV^(r,t2))(V^(r,ti)V'(r,ti)) 

(2.25) 

+ 4V.(a,(VV;(r, t2))^{r, t2)^{r, ti)V(r, ti))) , 

where o"^ is the Pauli matrix with entries aij = 1 — 5ij. The interaction Hamiltonian may now be 
obtained by integrating over all position space and over ti and t2 with < ti < t2- Note that the 
total divergences in (|2.24j) and (|2.25|) simply yield a boundary term on integrating with respect to 
r. We therefore have 

/roo i-t'z 
d-^r / dt2 / dti(V^(r,f2) 
Jo Jo (2.26) 

■ViP{r,t2)){^{r,h)iP{r,h)). 

We see immediately, by comparison with the free Hamiltonian TCq in (|2.8() . that the effect of this 
interaction is only to modify the diffusion constant. Thus if the trail has already visited a particular 
region of the lattice, it will experience a change in diffusion speed when it revisits this region. The 
nature of the modification depends solely upon the relative values of the bond and node interaction 
strengths, through the factor V' — V. The interaction derived above leads to a coupling proportional 
to Pi.(pi + q)- Therefore, referring back to the generalised interaction term displayed in (|2.1()|) . we 
see that in order to obtain the correct coupling for the kinetic self-avoiding trail on the Manhattan 
lattice, we must choose initially §2 = 0{V, V), and gi = g-i = 0{V'^, V'^, VV). This is in contrast 
to the kinetic self-avoiding trail considered by Peliti, for which gi > and (^2 = 53 = at t = 0. 
We see from (|2.26|) that, in the Manhattan case, 52 is proportional to V' — V, to first order. 
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2.2.3 Renormalisation group flows 

If 53 = initially, then to first order, ua = always. To simplify matters, we assume that 
this remains the case: this does not affect our qualitative conclusions. If we also choose ui = 
initially, then from Eq. 12.121 we find that the diffusion constant renormalisation factor Z = 1 at 
this point. This indicates that there is no singular contribution to the diffusion constant from the 
renormalisation at weak coupling, and therefore the dynamic exponent z = 2, with the possible 
appearance of corrections from the strong coupling behaviour; that is, time t scales as the square 
of the spatial extension r = |r(t)|. From the analysis presented earlier in this section, in d = 2 
dimensions the renormalisation group equations for ui and U2 in the M3 = plane are given by 

^ = ^ul + \u\ - 2uiU2 (2.27) 

-^ = -u\-ul + 3nin2 (2.28) 

Hence the only fixed point for e = is the trivial fixed point ui = U2 = 0. This corresponds to 
the case p = 0, for which the trails are necessarily non- interacting. The renormalisation group flow 
equations may be integrated to yield the general solution 

ui-U2 = A{ui + na)^/^ {2ui - ua)^/^ , (2.29) 

where if initially ui = and U2 = U2, then A = l/u2- Much may be learned by examining these 
flow trajectories. Fig. |31 illustrates schematically renormalisation group flow in the plane ^3 = 0. 
The attractors U2 = Mui may be identifled by dividing (|2.27l) by (|2.28j) and solving the resulting 
equation for M. There are three solutions: U2 = —ui, U2 = 2ui and U2 = ui. 

If asymptotically U2 = Mui for some M, then in this limit the renormalisation group equations 
give ui, ii2 ~ T~^. The flows in the direction of decreasing r in Fig. Olshould show the approach to 
asymptotic behaviour: see ^Hl- We observe that a trajectory starting at ui ~ and Ua > flows to 
the line U2 = —ui, whereas a trajectory beginning at ui ~ and U2 < flows to U2 = 2ui. In both 
cases, we obtain runaway with a flow to strong coupling. These two basins of attraction are divided 
by the separatrix U2 = ui, along which the flows also go to strong coupling. Although in all three 
cases the flows go out of the region of validity of the one-loop approximation, nevertheless their 
topology is consistent with what we might expect for the phase diagram of trails with attractive 
on-site interaction, depending on the relative strengths of the bond repulsion V and the node 
attraction V' . Recall that, to first order, U2 (x V — V'. For V ^ V' we would expect the walks to 
be in the universality class of trails, which is the same as that of ordinary self-avoiding walks pTlllHj . 
We may therefore interpret flows along the attr actor U2 = —ui as going towards a strong-coupling 
flxed point which corresponds to the universality class of SAWs, even though this point is beyond 
the range of applicability of our analysis. Similarly, for V ^ V we might expect the walks to be 
collapsed, so that we may interpret the attractor U2 = 2ui as flowing towards a (trivial) strong- 
coupling fixed point representing compact objects. Within this interpretation, then, flows along 
the separatrix go towards a strong-coupling flxed point corresponding to the collapse transition. 
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Figure 3: Schematic renormalisation group flow diagram in the plane U3 = 0. The arrows indicate 
the direction of decreasing r, i.e. increasing k. The three attractors are labelled (A) U2 = —ui, 
(B) U2 = 2ui and (C) U2 = ui. Possible trajectories are illustrated schematically by dashed lines. 
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that is the theta-point. (We note that for this simple phase diagram we need to assume that the 
initial value of ui, computed to second order in V and V', is slightly negative, otherwise there will 
be an intermediate phase in which the flows go to weak coupling, as for the Peliti model.) 

Although this analysis is based on a weak-coupling expansion, the topology of the consequent 
phase diagram should extend to strong coupling, unless there are further unexpected phase bound- 
aries. It does not, however, tell us the form of the phase boundary for larger values of V and 
V , and in particular in which phase the original lattice model lies. Thus we turn to numerical 
simulations in the next Section. 

Although the flows in (|2.27|2.28j) always go to strong coupling in our case, if we assume that 
they go to a fixed point describing a phase with a finite length scale (for example the radius of 
gyration in the collapsed phase), we may still use them to predict the form of the dependence of ^ 
on the initial, weak-coupling, parameters. This is because the largest contribution comes from the 
neighbourhood of the weak-coupling fixed point. In general, any such length scale satisfies 

e(nO) =eoexp dr^ , (2.30) 

where is of the order of the microscopic cut-off, r is the RG 'time', and r* is chosen so that 
^(u2(t*)) = 0(1). The integral here may be rewritten as 



dr , r^^^) du2 

A = / 2 2^ Q 

dit2 —u1 — ui + 6uiU2 



2 



Since initially ui = ©(u^), the integral is 0(1/^2) as U2 — > 0, so that ^ ~ exp(l/u2). This 
result would be academic were it not for the fact that the diffusion constant behaves as as 
p — > 0, and this places us squarely in the weak-coupling regime. Denoting the dimensionality of 
a quantity A by [A\, we find from 1)2. 8() that [^{T,t)ip{T,t)\ = k'^ and [D] = fik^'^, where /i is 
the Laplace-transformed temporal variable and k represents momentum. From (|2.26|) . we see that 
[52] k'^'^k'^ = fi'^k'^, i.e. [52] = ^J?k~'^~'^ . Hence 52 ~ D^k''. Therefore the renormalised, dimensionless 
coupling U2 ~ g2D~^K~''. Now since, from section [2.11 the bare diffusion constant D ~ , we 
have ~ p^. This leads to the leading order prediction ^ ~ e'^°°^* /P^: the correlation length 
increases exponentially with decreasing p^. This is entirely consistent with the only critical point 
of the lattice model being at p = 0. 



3 Simulations 
3.1 Preliminaries 

In this section, we use Monte Carlo simulations to study various characteristics of kinetic self- 
avoiding processes on the Manhattan lattice. The lattice bond length is chosen to be two units, so 
that a = 1 in the analysis above, and the lattice size is effectively infinite. The walker moves from 
one bond centre to an allowed neighbouring bond centre at each step, subject to the lattice bond 
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The vertical axis carries a logarithmic scale. 



directions and the condition that the path be self-avoiding; the number of steps taken is denoted 
by A^. Thus a path consisting of one step, for which = 1, connects two adjacent bond centres 
and has end-to-end length two units. Each simulation consists of a batch of 10^ walks, which are 
grown individually. When a walk reaches the specified maximum length or closes, that is, returns 
to its initial coordinates, it is terminated and a new walk is started; in this way, the walks sampled 
at a given length N should be statistically independent. Data points for different values of are 
generated from separate simulation runs. Errors on each data point are calculated as the variance of 
the quantity in question. Where curves have been fitted to the data, the fitting procedure employed 
an implementation of the nonlinear least-squares (NLLS) Marquardt-Levenberg algorithm. 

3.2 Survival probability 

Survival probability P{N) is defined as the proportion of simulated paths which have not yet closed 
after the A^th step. This is expected to be a j?-dependent quantity: at p = 1, P{N) = for A^ > 4, 
and at p = 0, P{N) = 1 for all A^ as the trails simply proceed straight along a given direction. The 
data plotted in Fig. EJsuggest that P(N) decays exponentially with increasing A^, on ap-dependent 
length scale. This corresponds to localisation for all non-zero p. Since all paths are extended at 
p = 0, this provides evidence in support of the existence of a trivial critical point, as discussed 
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A{p)-^ P{N/N*,p) exp{C{p)^/N/N*{p)), are plotted against the rescaled path length N/N* 
for selected values of The dashed line shows f{N/N*) = e-xp{—N/N*). 
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p 


N*{p) 


C{p) 


0.7 


34.0 


0.346 


0.6 


158 


1.15 


0.55 


338 


0.748 


0.5 


1320 


1.23 



Table 1: Values for coefficient C{p) and typical path length N*{p), tabulated for various j). Accuracy 
is to 3SF. 

earlier. In fact, we find that the simulated values for survival probability are consistent with a 
scaling form given by 

PiN/N*,p) ~ e-^/^*(^') e-^'^f^ ^N/N'{p) ^ (3^32) 

where N*{p) may be interpreted as the typical number of steps taken by a path before closure. 
Estimated values for C {p) and N* (p) , obtained by fitting the scaling form of 1)3.32(1 to the simulated 
data, are presented in Table ^ for certain values of p. Note that the typical length N* increases 
as p decreases; for this reason, it was found to be extremely difficult to obtain reliable estimates 
for N*{p) and C{p) for p < 0.5. Indeed, at p < 0.5 it is not possible to give reliable estimates 
for A^* since the estimated values are larger than the maximal value of used in the simulations. 
However, this is to be expected since as p — > 0, A^* should tend to infinity. Fig. ISlshows rescaled data 
A{p)~^ P{N/N*,p) ex.-p{C{p)^J N/N*{p) plotted against N/N*{p), where A{p) is the normalisation 
factor required for ()3.32|) . The data collapse on to the line f{N/N*) = exp{—N/N*), also plotted 
in Fig. [3 showing excellent agreement with the postulated scaling equation. Notice that since the 
typical length scale for closure of walks is much longer for smaller values of p, as predicted in the 
final part of section [2.21 data for these simulations is plotted over a smaller range of N/N*. The 
scaling behaviour of (|3.32|) may be justified by the following argument, which is motivated partially 
by evidence presented in section 13.41 for scaling of internal energy. 

The kinetic self-avoiding processes in which we are interested are non equilibrium problems. 
However, it is possible to use statistical mechanics to examine their structure by exploiting an 
equivalence with certain equilibrium objects. To demonstrate this equivalence, we employ argu- 
ments similar to those used by Bradley for kinetic growth walks on the Manhattan lattice jl9l I20j 
and by Owczarek and Prellberg for trails on hypercubic lattices \21\ 122) . Consider an equilibrium 
model whose configurations are kinetic self-avoiding trails laid down on the Manhattan lattice. For 
the moment, we look only at the case p = 1/2. The temperature is fixed at a particular value, as 
is usual with Monte-Carlo simulations. Now suppose that each doubly occupied node on a trail is 
assigned a weighting w. w = 1, then this weighting will simply yield ordinary trails, on which 
a node may be visited any number of times. However, if w is large then the path should be a 
highly compact object. For our kinetic self-avoiding trails, the appropriate weighting is w = 2: 
when the path arrives at a node for the first time, it may choose either of two possible exits, with 
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probability p = 1/2 of making each choice, but on revisiting the node, the path must leave by the 
one unoccupied exit with probability 1. Let 4>{N) be an open path configuration with length 
steps, on which n{(j)) nodes are doubly occupied. The probability of obtaining this configuration is 
P4i{N) = (1/2)^^"''''^^ However, as we have seen, there is also the possibility of generating a path 
which is a closed loop; the probability P{N) of obtaining an open configuration at length N is less 
than 1. In fact, 

HN) ^(N) (3.33) 

= 2-^ ^^(2), 

where 2^Ar(2) is the partition function for the equilibrium model, with the Boltzmann weight w = 2. 
Thus the set of kinetic self-avoiding growth trails at p = 1/2 is a sample from the Boltzmann dis- 
tribution of equilibrium trails with weighting w = 2. Further, the mean number of doubly occupied 
nodes is given by wdF/dw evaluated at tt; = 2, where F is the free energy of the equilibrium 
model. This relationship is particularly useful since the mean number of doubly occupied nodes 
may readily be computed from simulations. 

For a free random walk, Z]\f = 4^, where N is the number of steps taken, since the walker may 
leave a given node via any of the four associated bonds. Thus, in this case, F ~ A^log4. Note 
that since the bulk fractal dimension of a walk model is given by df = 1/z^, where the radius of 
gyration scales with exponent i^, here dj = 2. However, the free random walk does not have a well 
defined surface. In contrast, there are other types of walk for which a surface may be identified. 
For collapsed polymers, for example, which may be modelled by certain varieties of compact walk, 
a term appears in the partition function Z, and hence also in the free energy F which scales 
according to the surface area of the object [23^. For an object embedded in d dimensions, the 
average extension in any given spatial direction is given since v = 1/d. Therefore the 

surface area must scale as R'^^^ ~ ]\[{d-i)/d^ Thus in two dimensions, the surface area of a compact 
object scales like A^^/^. This leads to an expression of the free energy of the form 

F ~ iVFfc + N'^/'^Fs + . . . , (3.34) 

where the . . . indicates smaller correction terms, and F^ and Fg are respectively bulk and surface 
contributions to the free energy. Hence the partition function must scale as 

Zn{2) ~ exp(iVFi, + N'^I'^F, + ...). (3.35) 

From (|3.33|) . we may therefore deduce that if the kinetic self-avoiding trails behave as compact 
walks, the probability of finding an open configuration after steps should scale as 

P{N) ~ exp(iV log 2 + iVFfe + N^'^F, + ...). (3.36) 

This is consistent with the scaling form of 1)3.32(1 which was fitted to the simulated data, suggesting 
that the kinetic self- avoiding trails are compact objects. It follows from our postulation of a trivial 
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critical point that the scaling form ()3.36|) put forward for the case p = 1/2 should hold for all 
non-zero p. 

We note that our results indicate that kinetic growth trails on the Manhattan lattice are mapped 
to a temperature inside the collapsed phase of the equilibrium model. Previous kinetic growth 
models of walks and trails I2()| I21L |221 in two and higher dimensions have all mapped to 
temperatures in the extended phase or 0-point of the corresponding equilibrium model. This makes 
our kinetic growth model of interest from the point of view of polymer models. 

3.3 Radius of gyration and end-to-end distance 

In this section we consider scaling of the spatial extension of the kinetic self-avoiding trails. We 
begin by examining the mean square radius of gyration, {Rq), which is expected to scale as N'^'^, 
for some z^. For a non-interacting random flight, z/ = 1/2: {Rq) = No? (l — A^~^) /6, where a is the 
(constant) step length |2l]. It is also known that u = \/d for collapsed polymers in d dimensions. 
Fig. 121 shows simulated data for trails which are still open after N steps. On fitting a scaling 
form {R^) = A(jp) A^^'^^p) to each data set, we find that the estimated values for the exponent v 
lie in the range 0.47 ^ z/ ^ 0.52. Of course since the trails can never visit any bonds twice one 
can easily show that u > 1/2. We have also studied the mean square radius of gyration for walks 



17 



300 



A 250 



CD 
CD 

DC 

V 

Oi 

o 

CO 



200 



« 150 
■o 

T3 


o 100 

I 

T3 
CD 



CO 
DC 



50 







p=0.7 — ^ 
p=0.6 - X 

p=0.55 * 
p=0.5 

p=0.45 * 
p=0.4 « 
p=0.3 ^ ► 
p=0.2 ^ 







,06 



e © 



e e 



e e 



e e 



e © 



e e © «^ ' 



, © e © 



© © © 



,©' 




M « 



X X X ^ 



^ 5^ 5S ^ ^ 8« ' 



200 400 600 800 1000 1200 1400 1600 1800 2000 

Number of steps, N 

Figure 7: Scaling of RMS end-to-end distance with number of steps, for open trails, at various p. 
A line fitted to the data for p = 0.5 has equation {RD w 0.104A^ + 135iV^/^ - 728, showing hnear 
scaling with corrections 0(A^^/^). 



which close at the A^th step; the simulated data support a scaling form {Rq) ~ A, although the 
statistics obtained from these walks are not high. For example, for closed trails at p = 0.5, we have 
approximately {Rq) oc N^'^^. This evidence strongly suggests that for both open and closed trails, 
{Rq) ~ N'^'^ with 1/ = 1/2. This implies that the trails should be classified either as non-interacting 
random flights (polymer chains) or as compact walks (polymers in the collapsed phase) since both 
these have = 1/2 in two dimensions. However, evidence from the previous section suggests that 
the compact hypothesis is the most consistent one. 

We next consider the scaling of end-to-end distance Re for the trails. For the free random flight 
model, (i?g) = Aa^, where, as before, the step length is given by a. Thus for large A, {Rq) ~ {Rl)/6 
in this model. For compact walks, however, there should be corrections to this scaling form which 
result in suppression of the end-to-end distance relative to the non-interacting case. Fig. [3 shows 
simulated data for the root mean square end-to-end distance (i?g(A)) for trails at various values 
of p. Since this quantity does not make sense for paths which have closed, the average is taken 
only over those paths which are still open after the Ath step. We find that the data are consistent 
with a scaling form y/{Rl{N)) ~ N"^" (1 + BN-'^/'^ + CA^-i), where i/ = 1/2 as above. That is, we 
observe linear scaling with corrections of order A^/^ and smaller. The presence of the correction 
terms is consistent with behaviour expected from compact random walks, in contrast to the exact 
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linear scaling which would obtain if the walks were free random flights. From the result ^ ~ e^l^"^ 
in section 1221 we expect to be able to conclude that B^/N ~ e^^^^. However, this is difficult to 
verify. 

3.4 Internal energy 

In the equilibrium model described in section we may define the average internal energy U of 
a configuration of N steps by 

Q (3.37) 

since, up to constants, Fjv ~ logZ^r. Since Ff, and in 1)3. 34() are expected to be temperature 
dependent, taking the derivative of F with respect to temperature gives 

u ~ life + UsN~^/^ + ... (3.38) 

for the internal energy per step, U/N. Here, the corrections, which have carried through from the 
expression for free energy, are expected to be of the order of A^~^. The subscripts b and s denote 
bulk and surface terms, as previously. 

We now compare these results with simulations conducted on the Manhattan lattice. Fig. |H1 
shows data for internal energy per step, u{N), generated from populations of open trails of length N 
steps at p = 1/2. Application of a scaling form of the type given in (|3.38l) indicates that these kinetic 
trails exhibit behaviour consistent with the existence of a well-defined one-dimensional surface. We 
find the asymptotic behaviour u ~ u;, ~ 0.38, with Ug ~ 2.7, and as predicted, the corrections to the 
bulk and surface terms are of the order of A^^^. Since the bulk fractal dimension for these objects 
is l/i' = 2, we conclude that kinetic self- avoiding trails on the Manhattan lattice are compact, 
rather than extended, with a form similar to that of a liquid droplet. This conclusion is supported 
by the scaling forms already presented for survival probability, end-to-end distance and radius of 
gyration. Although the above analysis was carried out for the particular case p = 1/2, the same 
conclusions should hold for all trails away from the fixed point p = 0, since all trails generated at 
p > belong to the same phase. 

4 Conclusions and Discussion 

We have presented a model for quantum localisation in class C, on the Manhattan lattice. The 
general results of 7^ show that this may be mapped to a classical random walk problem, of trails 
with an attractive on-site interaction. The model depends on a single parameter p. A weak- 
coupling field-theoretic renormalisation group analysis appropriate for small p, numerical results 
for moderate values of p, and a rigorous argument for p > ^ , all lead to the conclusion that for 
all p > the walks almost always close after a finite number of steps, and correspond to compact 
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objects with a typical linear size which is a monotonically decreasing function of p. As p ^ 0, this 
size diverges hke exp (const. /p^). 

This picture leads to the conclusion that the quantum model on the Manhattan lattice is always 
in the localised phase, in contrast to the case of the L-lattice. This might have been expected, given 
the absence of edge states in the former case. We note that the weak-coupling RG equations for the 
classical walks have a critical dimension of two, as is the case for Anderson localisation. It would 
be interesting to compare the RG beta-function along the attractor U2 = 2ui with that obtained 
using sigma model methods This would involve relating our coupling U2 to the conductance, 
which is the physical coupling of the sigma model approach. 

The arguments above directly apply to the dynamic model of kinetic growth trails on the Man- 
hattan lattice. However, an equivalence exists to a particular temperature of a (static) equilibrium 
statistical mechanical model of polymers on that lattice, namely self-interacting trails. Hence, the 
results allow us to conclude that this temperature lies in the collapsed phase of the equilibrium 
model. This is the first time a kinetic growth walk/trail model has mapped into such a phase. 

There are several possibilities for extension of this work. The mapping established in j7j requires 
each node to have coordination number 4, with two incoming bonds and two outgoing bonds. This 
means that every node must either be of the type found on the L lattice, where a random walker 
is required to make a turn either to the left or to the right, or of the Manhattan type, where the 
walker may either make a turn or continue straight on. In this paper, we have considered only 
graphs consisting of nodes of one type; however one could also consider graphs consisting of an 
arrangement of nodes of both types. This work would be interesting in terms of both the quantum 
localisation behaviour and the nature of the corresponding classical walks. In jjj, Sp(4) bilayer 
systems were discussed, consisting of two coupled lattice layers each with Sp(2) symmetry. Work 
on Sp(4) models of this type is in progress this vein, one could examine mixed bilayer 

systems, in which the upper and lower layers have different lattice structures. For example, the 
localisation behaviour of a model with one L lattice layer and one Manhattan lattice layer might 
be expected to differ from that of a system consisting of two L lattice layers. The results presented 
in this paper and in [7] relate to models based in two dimensions. However, another possibility 
for further work would be to consider models whose underlying graphs are embedded in three or 
more dimensions. For example, in three dimensions, one could study a layered system, such as the 
three-dimensional U(l) model constructed by Chalker and Dohmen |26j . For systems in three or 
more dimensions, both insulating and metallic phases should exist; it is expected that each type of 
behaviour should hold over a range of values of p, with a transition from one to the other at some 
critical value p = Pc, where Pc G [0, 1]. It would be interesting to discover whether this is indeed 
the case. 
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